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Entanglement and violation of Bell inequalities are aspects of quantum nonlocality that have 
been often confused in the past. It is now known that this equivalence is only true for pure states. 
Even though almost all the studies of quantum correlations at quantum phase transitions only deals 
with entanglement, we here argue that Bell inequality can also reveal a general quantum phase 
transition. This is also shown for a particular case of two spin-^ in an infinity one dimensional chain 
described by the XXZ model. In this case, the Bell inequality is able to signal not only the first order 
phase transition, but also the infinity order Kosterlitz-Thouless quantum phase transition, which 
can not revealed neither by the energy of the system nor by the bipartite entanglement. We also 
show that, although the nearest-neighbors spins are entangled they, unexpectedly, never violate the 
Bell inequality. This indicates that the type of entanglement which is relevant for quantum phase 
transition is not trivial, i.e., can not be revealed by the Bell inequality. 



I. INTRODUCTION 

The concept of entanglement has evolved since it has 
first appeared in 1935, as a "spooky action at a distance", 
i.e., the possibility of one system to influence another in- 
stantaneously, even at large distances Only in 1964, 
John Bell moved the concept from the philosophical de- 
bates to the laboratory, using inequalities - now known 
as Bell inequalities [2,]. Entanglement was then a prop- 
erty of states which could not be described by realistic 
local theories. Perhaps the more drastic evolution of the 
concept came in 1989, with the work of R. F. Werner Q. 
Werner defined entangled states as the ones which could 
not be created using local operations on the systems with 
classical communication between the parts; these opera- 
tions may create only classical correlations. However, 
Werner surprisingly found mixed states which were en- 
tangled, according to his definition, but did not violate 
the Bell inequality. Entanglement and non-locality be- 
came distinct, but related, concepts. 

Even though entanglement has all these facets, the 
word is usually used to designate quantum correlations 
that are "stronger" than classical ones. On the other 
hand, phase transitions {e.g., magnetic phase) are gener- 
ally related to an emerging order with origin in long range 
correlations. Quantum phase transitions (QPT) are the 
ones occurring at zero temperature, where the systems is 
in a pure quantum state. Therefore, we may expect these 
long range correlations to come from entanglement. 

Studies trying to link entanglement and QPT first ap- 
peared in 2002 for specific models 0, [^. Later, it was 
shown that, in principle, entanglement should inherit 
the non-analytical behavior of the ground state energy 
at any QPT; since both originate in the ground state 
0i0|- However, accidental non-analytical behavior may 
also appear, or disappear, due to the maximization pro- 
cess involved in entanglement definitions (sj. Nowadays, 
there are many studies concerning the link between en- 
taglement and QPT. Their results can be favorable or 
unfavorable to it, depending of the specific model and 



the measurement of entanglement used (for a review see 
0). Nonetheless, to find such a link, all these studies 
used the Werner definition of entanglement and did not 
analyze non-local properties of the states. 

In this work, we complement the study of the relation 
between entanglement and QPT, analyzing the non-local 
aspects of quantum phase transition using Bell inequali- 
ties. The objective is to find: i) if entangled states, that 
may play a role in a quantum phase transition, are also 
non-local; ii) if the Bell Inequality is able to signal QPT. 
We argue that, in principle, the Bell Inequality should 
also signal quantum phase transitions. A particular case 
of two spins in an infinite one-dimensional spin-i XXZ 
Heisenberg chain was then studied. The study of non- 
local aspects of QPT is very recent, with only one analysis 
of a particular case of the one dimensional XY Heisenberg 
chain with a transverse magnetic field published[10][35]. 

In order to be self contained, we organize the article 
as follows: In Sec. II we give a brief introduction about 
entanglement and non-locality through Bell Inequality, 
highlighting its differences. Sec. Ill describes the XXZ 
model and its relevant quantities, in particular how to 
compute its correlation functions. The mathematical ex- 
pression for the entanglement and Bell Inequalities in the 
specific case of the XXZ model are given in Sec. IV. The 
main results and a detailed discussion of their implica- 
tions are in Sec. V and we conclude in Sec. VI. 



II. ENTANGLEMENT AND NON-LOCALITY 

Bell inequalities investigate correlations between two 
parties who share a quantum state. In the case of two 
spins 1/2 particles, each side, A and B, choose one direc- 
tion, di and bj, and measures, simultaneously and inde- 
pendently, the observables Ai —Oi.a and Bj =bj.a; with 
a — (fj^, cr^, (7^). If each part picks two directions. Bell 
inequalities state that any theory which is both realistic 
and local gives [s^ 
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III. SPIN-i XXZ MODEL 



\{Bchsh)\ =- |(Ai ® Bi) + (Ai (E) B2) + (A2 ® Bi) 
-(A2®B2)|<2. (1) 

The value obtained for (Bchsh) depends both on the 
state and on the directions chosen. So, for a given state 
one should maximize over all directions defining a Bell 
measurement as: 



B = max I {Bchsh} 



(2) 



The states for which B > 2 are the ones which cannot 
be described by a realistic local theory and were then 
associated with non-classical correlation or non-locality 
[syj . Clearly, states which can be written in a product 
or separable form, are local. These states are 

known as separable or not entangled. 

In 1989, Werner [3] proposed a definition of separabil- 
ity for mixed states from an operational point of view: 
separable states are the ones which can be produced by 
local operations and classical communication (LOCC). 
They can be written as 



k 



(3) 



with < pfc < 1 and X^fcPfc ~ 1- Separable mixed states 
can be correlated, but this correlation is a classical one 
that comes from the probabilities pk. Surprisingly, in the 
same work, Werner constructed states which are entan- 
gled by his definition, but nonetheless do not violate the 
Bell inequality. Entanglement, although necessary, is not 
a sufficient condition for a mixed state to be nonlocal; 
Bell inequality does not capture all "kinds of entangle- 
ment". 

For pure bipartite states, only in 1991, Gisin proved 
that entanglement is a necessary and sufficient condition 
for a state to be non-local [Uj. Until nowadays, there 
are many attempts to draw a unifying picture between 
non-locality and entanglement by generalizing the Bell 
inequality. Popescu, in 1995, for example, showed that 
an entangled state not violating Bell, could be made to 
violate it after some processing with LOCC [l3|- Another 
possibility, is to allow each part to make a measurement 
in a third direction. This originates a new inequality 
that is violated by states which do not violated CHSH 
inequality |13l|. F or a introduction and review on similar 
attempts see |14| . There are also many other definitions 
of non-locality, and all of them show some anomaly with 
entanglement fl5|. For example, while for two qubits 
pure states the degree of violation of the Bell inequality 
is a increasing function of the entanglement of the state, 
this is not true for two qudits. Finally, there are also 
other definitions of quantum correlation, as the discord 
which has also been studied at quantum phase transitions 



We now describe the physics of the model we intend 
to use to study the relationship between entanglement 
and non- locality: an one dimensional spin-^ chain where 
the spins interact through anisotropic Heisenberg inter- 
action. It is known as XXZ model and for N particles its 
Hamiltonian is given by 



N 



qy qy 



(4) 



where S'lj = cr"/2 [u = x,y,z), cr^ are the Pauli spin-i 
operators on site j, A is the anisotropy parameter, and 
'^j+N — '^j • The XXZ model cannot be diagonalized, but 
its energy spectrum can be obtained by Bethe ansatz. 
The Hamiltonian has two symmetries: i) a discrete par- 
ity Z2 symmetry over the plane xy: ct^ — >■ — cr^, and ii) a 
continuous C/(l) symmetry that rotates the spins in the 
xz plane by any angle 6 . The Z2 symmetry implies that 
(erf) = and {afa]) = {a^ap = 0, while the C/(l) sym- 
metry implies that (af) = (erf ^ — 



and 



0, (afap = (afaj) 



. However, in the thermodynamical limit 



a quantum phase transition occur and breaks the discrete 
Z2 symmetry; a continuous symmetry may not be broken 
at one-dimension even at zero temperature. The model 
has three phases: 

i) A < — 1: the system is in a ferromagnetic phase with 

all the spins pointing in the same direction. There 
is a first order quantum phase transitions (IQPT) 
at the critical point A = — 1 . 

ii) — 1 < A < 1: the system is in a gapless phase, where 

the correlation decays polynomially. 

iii) A > 1: the system is in the anti- ferromagnetic phase. 

The critical point at A = 1 is of infinite order 
or Kosterlitz-Thouless quantum phase transitions 
(KT-QPT). 

The Bethe ansatz solution gives the ground state energy 
MM: 



eo(A) 



A 
' 4 ' 



A I sin niy r° 
4 ' 2tt J_ 

3 -In 2, 



dx 



sinh X sinh ux 



A<-1, 
,-l< A < 1 , 
A = l, 

(5) 

For A > 1 we just need to change 



where A = costt:/. 
— in (6). 

The correlation functions, on the other hand, cannot 
be easily obtained from the Bethe ansatz solution. They 
are given in terms of very complicated multiple integrals, 
and solving them is a non trivial mathematical problem. 
However, for nearest neighbors we can obtain the corre- 
lation from eo(A): 



+1/ 



. deojA) 
dA ' 



(6) 



3 



(<^r<+i) = (^r^f+i) = ^(4eo(A) - A(afaf+i)). (7) 

For spins further apart, progress has been slow, but there 
are already some expressions available up to third neigh- 
bors [2Qi] • We will not show them here, since they are too 
lengthy J^. 

We are interested in the entanglement and non-local 
properties of two spins in the chain. These can be de- 
termined by the reduced density matrix of the two spins, 
which can be obtained from the magnetizations and cor- 
relations. Actually, the density matrix of any two spins- 
ican be expressed as 



Pi{i+r) — 



I (g) I -I- p.o- (g) I -I- I (g) q.o- 4- ^ f^o-r ® 



(8) 

with r being the distance between the sites, p = (ct (X) I), 
q = (I (g) (t), and t^f" = {(Tfaf_^_^). Due to the Hamilto- 
nian symmetries, only {i^^, t^^} do not vanish and 
= t^'^. Translation invariance also implies that the 
correlation functions for two-sites depend only on the 
distance between the sites (r), being independent of i. 
So for the XXZ model we have 



Pr 



+ 



V 





1-tf. 







2i^^ 

1 - tf 

1 + 



(9) 



Remark that in the ordered phases, ferro and anti- 
ferromagnetic, due to the spontaneous symmetry break- 
ing (SSB), a finite value of the magnetization will emerge: 
(erf) = TO. This should be taken into account, complicat- 
ing a little bit the expression for p. However, in many 
studies this is not done, and the simpler expression, pre- 
serving the symmetries, is used. In this case, the ground 
state used is not the one which breaks the symmetry 
and is stable against external perturbations, but a su- 
perposition of the two degenerated ground states with 
opposite magnetizations. This could be justified for the 
study of very small systems where the SSB does not hap- 
pen, but not for real macroscopic systems. It is true, 
however, that all the thermodynamical quantities, which 
comes from the energy eigenvalues, are not affected by 
this choice. Nonetheless, entanglement properties may 
be affected by this choice, since it also depends on the 
eigenvectors j2T|. [2^. Last, besides the symmetric super- 
position of the two ground states, one could also think 
of an equiproportional mixture of the two, that would 
emerge if one takes the limit of temperature going to 
zero. 



IV. BELL MEASUREMENT AND 
CONCURRENCE 

In general, checking a Bell inequality violation is a 
high-dimensional variational problem. We have to op- 
timize over the various possible measurements settings 



that each observer performs (see Eq. One sim- 

plification is that the value of Bchsh and B does not 
depend on local properties, the magnetizations p and q, 
but only on the correlations i"". For two spin-i particles 
a closed analytic expression for the maximum violation of 
the Bell inequality was obtained by the Horodeckis (23| 
and is given in in terms of the matrix U — T^T by 



(10) 



with u and u' the two largest eigenvalues of U and T the 
3x3 matrix build from the correlations t^"". As u and 
u' are combinations of the correlation functions, we can 
rewrite it in terms of the correlations. 

The expression for B in terms of the correlations is very 
cumbersome for general states. But, because of the sym- 
metries of the XXZ model, see Eq. ©, it becomes quite 
simple. In fact, the matrix T is already diagonal, with 
^xx ^ ^yy g^^^ ^zz j|.g diagonal elements and t^^ — t^^ . 
Therefore, the eigenvalues of U are (i^^)^ and (t^^)'^ and 
we only have to check which is greater. The expression 
becomes 



Br = 2 max{y^2(afaf+,)2, ^{afaf^^^ + («+,)2}. 

(11) 

To have non-locality, we need either 

• i) I («+.)! > |«af+,)| and I («+,)! > 



or 



ii)l(f^f<+r)l 



< 



\2 



'i+r/ 



'i+r/ 



2 + r y 



> 1. 



To compare the non-locality with entanglement between 
the two spins, we will use the concurrence as our en- 
tanglement measure; see (24| for a review on entan- 
glement measures. It is a well defined entanglement 
measure and is given in terms of the spectrum of pp', 
where p' ^ ® (7^ p*ay (E) is the time- reversed den- 
sity matrix. Let A; be the eigenvalues of pp' so that 
Ai > A2 > A3 > A4. Then the concurrence is 



(12) 



C ~ max{0, v Ai 



A4}. 



The concurrence can also be written in terms of the 
correlation functions, but in general the expression is also 
very cumbersome, since it comes from solving a polyno- 
mial of fourth degree. As in the case of the Bell in- 
equality, due to the symmetries, for the XXZ model the 
expression is simple and given by (25j 

r 2|«af^,)|-(l + «a,f^,)) ] 
Cr = max <^ 0, S . (13) 

As mentioned before, the Bell inequality is not directly 
affected by local properties, the magnetizations, but only 
depends on the correlations, given by the matrix T. The 
concurrence, however, may be affected, since the local 
magnetizations may appear in its expression. Eq. (|13p 
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is only valid when the magnetizations and off-diagonal 
correlations are null {{<jf) — and {ufu^) — (cfcrj) = 
(crf(T|) = 0), and no simple expression can be obtained 
in terms of the correlation for the general case. This 
may have consequences to the effect of the spontaneous 
symmetry breaking on these properties. In the case of the 
XXZ model, the SSB manifests itself with a finite value 
of (af), which could in principle change the value of the 
concurrence, but not of the Bell inequality. However, as 
showed in |25|, it turns out that for the XXZ model the 
concurrence does not change when the SSB is taken into 
account. 




V. RESULTS 

Before showing the results of non-locality and compare 
with the entanglement, let us look at the ground state 
energy per site eg, as shown in Fig. [T] We can see the 
non- analytical behavior of at the IQPT A = — 1 and 
the nonexistence of non-analyticities at the KT-QPT at 
A = 1. We also show the correlation functions for first, 
second and third neighbors in figs. [5] and [31 A non- 
analytical behavior can also be seen at A = —1, and at 
least for first neighbors, it is directly related to the one 
found in cq. 




Figure 2: Correlation function in the x direction, (afaf^,.), 
for first, second and third neighbors of the XXZ model as a 
function of the anisotropy A. It can be seen that the correla- 
tion is able to signal the first order quantum phase transition 
at A = — f, but not the infinite order quantum phase transi- 
tion at A = f. 
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Figure 3: Correlation function in the z direction, {aiolj^^), for 
first, second and third neighbors of the XXZ model as a func- 
tion of the anisotropy A. It can be seen that the correlation 
is able to signal the first order quantum phase transition at 
A = — 1, but not the infinite order quantum phase transition 
at A = 1. 



Figure I: Ground state energy per site of the XXZ model as a 
function of the anisotropy A. There is a first order quantum 
phase transition at A = — 1 from a ferromagnetic to a gapless 
phase. At A = 1 there is an infinite order Kosterfitz-Thouless 
phase transition from the gapless to the anti-ferromagnetic 
phase, which is not indicated by cq. 

Last, we review the results of the bipartite entangle- 
ment between two spins, first obtained in [26]. In Fig. 2] 
we observe that the entanglement of first neighbors sud- 
denly appears as they enter the gapless phase, achieve 
the maximum at the KT-QPT, and slowly decrease in 
anti-ferromagnetic phase. Second and third neighbors 
are only slight entangled and in a small region on the 
right of the IQPT; the curves can barely be seen in Fig. 
31 The first order phase transition is indicated by the 
non-analytical behavior of the concurrence of all the first 



three neighbors: its first derivative diverges (see Fig. [SI 
for the case of first neighbors) 

We now present our results. First we argue that the 
Bell Inequality may be as good as entanglement to sig- 
nal the QPT. The argument is equivalent and follows 
directly from the one used in Q for bipartite entangle- 
ment and in [27] for multipartite: for Hamiltonians with 
only nearest-neighbors interactions not only the entan- 
glement between nearest-neighbors spins, but also cq can 
be obtained from the reduced density matrix of the two 
spins; they are functions of Pi^i+i- So any non-analytical 
behavior in eo = f{pi,i+i) will also appear in the entan- 
glement E = g{pi,i+i)i unless the the specific form of the 
functions / and g cause the non-analitycal behavior or 
originate accidental cancellations (see [8] for example). 
Note that while / is a linear function of Pi,i+i, g involves 
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the absolute value of and maximizations or mini- 

mizations procedures. 

But such analyze is also trivially true for any property 
of the two spins, since they are completely determined 
by Pi^i+i, including the Bell measure. Therefore B is, in 
principle, as good as any entanglement measure to signal 
a QPT. But as in the case of entanglement accidental can- 
cellations may also happen due to the complicated form 
B depends on Pi^i+i'. it also involves the absolute value 
of pis+i and a maximization. Note that signaling the 
QPT is not the only interest when studying the behavior 
of entanglement or non-locality at QPT. One could also 
learn about the nature of the correlations involved in the 
QPT, as we will argue in the following when looking at 
the specific case of the XXZ model. 

Given the general argument, we move to our results for 
the case of the XXZ model. The value of the Bell mea- 
surement, using Eq. is shown in the upper curves 
of Fig. Hlfor first (solid blue curve), second (dashed red 
curve) and third (doted yellow curve) neighbors. We ob- 
serve that none of them violate the Bell inequality, de- 
spite being entangled. But lets first concentrate on the 
analytical properties of B. As it may be expected from 
the arguments above, it is able to indicate the QPTs. 
Furthermore, not only the first order QPT is revealed by 
B, but also the Kosterlitz-Thouless QPT, which is neither 
indicated by the energy nor by the bipartite entanglement 
given in terms of the concurrence. This can also be seen 
in Fig. [6] where we plot the derivatives of the Bell mea- 
surement Bi and the concurrence Ci for nearest neigh- 
bors. It shows that while the derivative of both diverges 
at the first order quantum phase transition (A = — 1), 
only the derivative of the Bell measurement signal the 
infinite order quantum phase transition (A = 1). 

It is interesting to note that the success of B in re- 
vealing the KT-QPT comes from the maximization pro- 
cess, which is usually associated with negative results 
in the ability of entanglement to signal a QPT. More 
specifically, analyzing Eq. (|lip . we find out that while 
in the ferro and anti-ferromagnetic phase Ktrfcrf^^)! < 

|(afaf+,,)| and Br = 2^(«+,)2 + (afaf_^,)2, in the 
gapless phase |(crfcrf+r)l > \{(^t<^i+r)\ Br = 

2^2{a^ afj^^)"^ . On the other hand, the non- analytical 
behavior of B at the first order QPT does originate in 
Pi^i+r or in eo. 

Just to contrast, we should mention that the non- 
analytical behavior of the concurrence does not come 
from the non-analyticities in eo, but from the maximiza- 
tion involved in the concurrence. Note also that at the 
KT-QPT the first neighbor concurrence is maximum, 
while all the next neighbors are null almost in the whole 
phase. However, even though one can show that, in gen- 
eral, the concurrence may inherit the non-analytical be- 
havior of eo , it is not possible to show that it should be 
maximum at the critical point. Therefore, such a behav- 
ior, should not be an indication of a QPT and may be 
specific to the XXZ model; in fact it has been interpreted 



as a consequence of the specific symmetries of the model 
(see Sec. IV-A-2 of [9]). It would be interest to know if 
the behavior of B at the KT-QPT is also only specific to 
this model. 

Last, let us comment on the non- violation of the Bell 
inequality for nearest-neighbors, our third result. It is 
unexpectedly, since these spins are entangled, and not 
only slightly. They have a kind of hidden entanglement, 
similar to the one present in the Werner states, that is not 
revealed by the Bell inequality, being in this sense local. 
That entanglement may be relevant for quantum phase 
transitions has been suggested by many works [9]. Our 
results complements such works, indicating that the kind 
of entanglement that is relevant for QPT is not trivial: 
it cannot be revealed by the Bell inequality. Note that, 
one may argue that this should be expected since we are 
studying two particles in a huge chain in a pure highly 
entangled state. Thus, these two particles are in very 
mixed states close to the identity and should not violate 
the Bell inequality. But this argument also applies to 
entanglement: the two particles should not be entangled. 
In our case they are, and not slightly. In fact, despite 
even slight entangled states may violate Bell, the num- 
ber of states violating Bell and the degree of violation 
increase with the amount of the entanglement [l^. We 
also observe that second and third neighbors do not vi- 
olate Bell, but these are only weakly entangled in a tiny 
region of the phase diagram. Note that, even tough our 
states have a behavior similar to Werner states, they are 
not Werner states. Only at A = 1 can our state be writ- 
ten in the form of a Werner state. Thus, these are a new 
class of entangled states, which do not violate the Bell 
Inequality. 
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Figure 4: Bell measurement (three upper curves) and concur- 
rence (three lower curves) for first, second and third neighbors 
of the XXZ model as a function of the anisotropy A. It can 
be seen that the Bell inequality is never violated, even when 
the two spins are entangled. While both the concurrence and 
the Bell measurement are able to signal the first order quan- 
tum phase transition at A = — 1, only the Bell measurement 
reveals the infinite order QPT at A = 1. 
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Our results can be visualized geometrically in Fig. 
[5l where we have the parameter space of density ma- 
trices of two spin-i particles with the symmetries of 
the XXZ Hamiltonian. Remember pi j only depends on 
{(j'f(j'j) and (afaj); Eq. The greater triangle defines 
the physical states and the inner diamond the separable 
states (SEP). The two smaller triangles, inside and on 
the lower vertices of the greater, are the entangled states 
(ENT). The two smallest regions, also on the lower ver- 
tices of the greatest triangle, are the states violating Bell 
inequality (NL) . The parallel lines are contour lines of the 
concurrence. We can see that it increases in the direction 
of the lower vertices of the greatest triangle. The contour 
lines of the Bell inequality are parallel to the lines defin- 
ing the border of the region NL, but not show. We also 
plotted the "trajectory" of pi,i+r for r = 1, 2 and 3 as the 
parameter A is varied. For all the curves we marked fours 
points that corresponds to A = — 1 (full circle), -0.999 
(square), (triangle) and 1 (diamond). The whole ferro- 
magnetic phase (A < — 1) is represented by the full circle. 
As we move from the ferromagnetic phase in the direc- 
tion of the gapless (— 1 < A < 1) there is a jump at the 
point A = — 1. The state of Pi^i+i (solid blue curve) and 
/3i^i+3 (dotted yellow) jumps through the left side, enter- 
ing the entangled region, while Pi,i+2 (red dashed curve) 
jump to the entangled region through the right side. We 
can also observe that and pi,i+3 barely enter the 

entangled region and go back to the separable one, while 
Pi^i+i goes deeper in the entangled region and stays there 
in the whole gapless phase. The maximum value of the 
concurrence for pi.i+i: at the critical point A = 1, can be 
seen as the point where the trajectory tangencies a con- 
tour line of the concurrence. At this point the trajectory 
also touches a minimum of the Bell measurement, which 
happens due to the optimization procedure involved in 
its definition. 



VI. CONCLUSION 

We have studied non-locality using Bell inequality at 
quantum phase transitions. We first showed that the 
Bell Inequality is, in principle, as good as any entan- 
glement measurement to signal a quantum phase tran- 
sistion. Then, we studied the particular case of two spins 
in an one-dimensional spin-^ XXZ infinite chain at zero 
temperature and found that: i) the Bell inequality is not 
only able to reveal the first order quantum phase tran- 
sition, but also the Kosterlitz-Thouless quantum phase 
transition. This can not be indicated by the energy of the 
system nor by the bipartite entanglement given in terms 
of the concurrence, ii) the two spins, despite being entan- 
gled, never violate the Bell inequalities: their correlations 
can be then described by a realistic local theory. 

Our results suggest that the Bell inequality may con- 
tribute to the study of quantum phase transitions, since 
they are able to signal quantum phase transitions - and 
even the Kosterlitz-Thouless quantum phase transition. 




Figure 5: Parameter space of the reduced density matrix of 
two spins with the symmetries of the XXZ Hamiltonian. The 
greater triangle defines the physical states. The two smaller 
triangles, inside and on the lower vertices of the greater, are 
the entangled states (ENT). The two smallest regions, also 
on the lower vertices of the greatest triangle, are the states 
violating Bell inequality (NL). The curves are the "trajectory" 
of the ground state of the XXZ Hamiltonian as the parameter 
A is varied for pi,i+i (solid blue curve), pi,i+2 (dashed red 
curve) and pi,i+z (dotted yellow curve). For all the curves we 
marked fours points that corresponds to A = — 1 (full circle), 
-0.999 (square), (triangle) and 1 (diamond). 
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Figure 6: Derivative of the Bell measurement Bi and the 
concurrence Ci in relation to A for nearest neighbors. 



in our specific model. Furthermore, they may increase 
our understanding of the nature of the quantum correla- 
tions involved in quantum phase transitions. Again, in 
our specific model, it was shown that the type of entan- 
glement which is relevant at the quantum phase transi- 
tion is non-trivial, since it cannot be revealed by the Bell 
inequality. Therefore, Bell inequality constitute a dif- 
ferent resource, other than entanglement, and could be 
an alternative and complementary way to characterize 
a quantum phase transition, maybe even infinite order 
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ones. 

Our work gives information about the type of quan- 
tum correlation involved at the quantum phase transition 
of the XXZ model. To better characterize these quan- 
tum correlations, it would be interesting to study other 
models and more general Bell inequalities involving more 
measurement directions. The study of multipartite Bell 
inequalities is also of great interest, despite its technical 
difficulties. First, obtaining correlations among three or 
more spins in condensed matter systems suffering quan- 
tum phase transitions is not trivial. Second, there is no 
unique multipartite inequality, but many to be studied. 
However, a work on finite chain up to six spins was re- 
cently published [2^. 
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